Entanglement Witnesses in Spin Models 
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We construct entanglement witnesses using fundamental quantum operators of spin models which 
contain two-particle interactions and have a certain symmetry. By choosing the Hamiltonian as 
such an operator, our method can be used for detecting entanglement by energy measurement. We 
apply this method to the Heisenberg model in a cubic lattice with a magnetic field, the XY model 
and other familiar spin systems. Our method provides a temperature bound for separable states 
for systems in thermal equilibrium. We also study the Bose-Hubbard model and relate its energy 
minimum for separable states to the minimum obtained from the Gutzwiller ansatz. 

PACS numbers: 03.65.Ud, 03.67.-a, 05.50.+q 



I. INTRODUCTION 

Entanglement lies at the heart of quantum mechan- 
ics and plays also an important role in the novel field 
of Quantum Information Theory (QIT, 0|). While for 
pure quantum states it is equivalent to correlations, for 
mixed states the two notions differ. In this general case, 
a quantum state is entangled if its density matrix can- 
not be written as a convex sum of product states. Based 
on this definition, several sufficient conditions for entan- 
glement have been developed U|. In special cases, e.g. 
for 2 X 2 (two-qubit) and 2x3 bipartite systems lj| and 
for multi-mode Gaussian states Q even necessary and 
sufficient conditions are known. 

However, in an experimental situation usually only lim- 
ited information about the quantum state is available. 
Only those approaches for entanglement detection can be 
applied which require the measurement of not too many 
observables. One of such approaches is using entangle- 
ment witnesses. They are observables which have a pos- 
itive expectation value or one that is zero for all separa- 
ble states. Thus a negative expectation value signals the 
presence of entanglement. The theory of entanglement 
witnesses has recently been rapidly developing [J| . It has 
been shown how to generate entanglement witnesses that 
detect states close to a given one, even if it is mixed or 
a bound entangled state [^j. It is also known how to 
optimize a witness operator in order to detect the most 
entangled states Q- 

Beside constructing entanglement witnesses, it is also 
important to find a way to measure them. For example, 
they can easily be measured by decomposing them into a 
sum of locally measurable terms . In this paper we fol- 
low a different route. We will construct witness operators 
of the form 



Wc 



O 



inf [(*|0|*)]. 



(1) 



where S is the set of separable states, "inf" denotes in- 
fimum, and O is a fundamental quantum operator of a 
spin system which is easy to measure. In the general case 
inf^g5(^|0|^) is difficult, if not impossible, to compute 



[g. Thus we will concentrate on operators that contain 
only two-particle interactions and have certain symme- 
tries. We derive a general method to find bounds for the 
expectation value of such operators for separable states. 
This method will be applied to spin lattices. We will also 
consider models with a different topology. 

If observable O is taken to be the Hamiltonian then our 
method can be used for detecting entanglement by energy 
measurement [^. While our approach does not require 
that the system is in thermal equilibrium, it can readily 
be used to detect entanglement for a range of well-known 
systems in this case. The energy bound for separable 
states correspond to a temperature bound. Below this 
temperature the thermal state is necessarily entangled. 
Numerical calculations have been carried out for some 
familiar spin models. They show that for the parameter 
range in which substantial entanglement is present in the 
thermal ground state, our method detects the state as 
entangled. Thus our work contributes to recent efforts 
connecting QIT and the statistical physics of spin models 

iia. 



II. 



ENERGY BOUND FOR SEPARABLE SATES 



We consider a general observable O on a spin lat- 
tice defined in terms of the Pauli spin operators tr^'') = 



{^. 



(k) (k) (k) 



,cr' 



V 



(Jz 



Oas 



0:=0 



{^«} 



k=l 



(2) 



where O is some multi- variable function 'Tl'|. We will 
discuss how to find the minimum expectation value of 
such an operator for separable states of the form 



P 



5]piPp^®Pp^®-®P; 



(N) 



(3) 



The minimum of (O) for pure product states is ob- 
tained by replacing the Pauli spin matrices by real vari- 

(k) . WTTl 

ables s^ , , in Eq. |(2|| and minimizing it with the con- 
straint that s^*^^ are unit vectors [IJ. The minimum ob- 
tained this way is clearly valid also for mixed separable 



states since the set of separable states is convex 



^scp 



:= inf (^lOI*) = inf O 



{^''ru=. 



(4) 



In the most general case many-variable minimization 
is needed for obtaining Ogcp- In some cases, to which 
many of the most studied lattice Hamiltonians belong, 
it is possible to find a simple recipe for computing the 
minimum of O. 

(i) Let us consider an operator O which is the sum of 
two-body interactions. It can be described by a lattice 
or a graph. The vertices V := {1,2, ...,7V} correspond 
to spins and the edges between two vertices indicate the 
presence of interaction. 

(ii) Let us assume that this lattice can be partitioned 
into sublattices in such a way that interacting spins cor- 
respond to different sublattices. Fig. Q] shows lattices 
of some common one- and two-dimensional spin models. 
The different symbols at the vertices indicate a possible 
partitioning into sublattices with the above property. For 
simplicity, next we will consider the case with only two 
disjoint sublattices, A and B, and assume that O can be 
written in the form 



-• — ■ — •- 



^>__li_H^ 



(a) 



(b) 





(c) 



(d) 



FIG. 1: Some of the most often considered lattice models: (a) 
Chain, (b) two-dimensional cubic lattice, (c) hexagonal lattice 
and (b) triangular lattice. Different symbols at the vertices 
indicate a possible partitioning into sublattices. 
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and 






(5) 



denotes 



where / is some two-spin function, 
spins of sublattice A/B. 

If conditions (i) and (ii) are met, then it is enough to 
find spins s and s^ corresponding to the minimum of 
/(s'^, s^). Then setting all the spins in sublattice A to s^ 
and in sublattice B to s^, respectively, gives a solution 
which minimizes O. 



III. EXAMPLES 



A. Heisenberg lattice 

Let us consider an anti-ferromagnetic Heisenberg 
Hamiltonian with periodic boundary conditions on a d- 
dimensional cubic lattice 



Hh-.^Y.^- 



{k)^{i) 



.^(fe)^(O 
"y V 



-a^al') 



Ba^J'l (7) 



{k,i) 



The strength of the exchange interaction is set to be J = 
1, B is the magnetic field, and {k,l) denotes spin pairs 
connected by an interaction. The expectation value of 
Eq. 10 for separable states is bounded from below 



In the following we will use the Hamiltonian H for con- 
structing entanglement witnesses. The energy minimum 
for separable states is the same as the ground state of the 
corresponding classical spin model. Our method detects 
entanglement if 



AE 



(H) 



Escp < 0. 



(6) 



If AE < then | A£^| characterizes the state of the system 
from the point of view of the robustness of entanglement. 
It is a lower bound on the energy that the system must 
receive to become separable. 

We will use the previous results to detect entan- 
glement in thermal states of spin models. In ther- 
mal equilibrium the state of the system is given px — 
exp{-H/kBT)/Tr[exp{-H/kBT)] where T is the tem- 
perature and kB is the Boltzmann constant. For sim- 
plicity we will set fcs = 1. Using Eq. ^ a temperature 
bound, Te, can be found such that when T < Te then 
the system is detected as entangled. 



{Hh) > -E_ff,scp : — 



-dN{B^I 

-dN{\B\ 



1) 



if |B| 
if IBI 



<4, 
>4, 



(8) 



where N is the total number of spins. This bound was 
obtained using two sublattices, minimizing the expres- 
sion /ff (s-^, s^) := s^s^ -f B{s^ + sf )/2. Based on this 
EH,s.p^dNini[fH] M 

Let us now consider a one-dimensional spin-1/2 
Heisenberg chain of even number of particles. Ii B — Q 
then Eq. |(SJ corresponds to 



{k,i) 



')>-!. 



(9) 



which is simply a necessary condition for separability in 
terms of nearest- neighbor correlations. In the large N 
limit, the energy minimum for entangled states can be 
obtained as £;„„ = -4iV(ln2-l/4) « -l.llN HI. The 
energy gap between the minimum for separable states 
and the ground state energy of Hh is thus AEgap ~ 



Ep{B,T) 




ground state energy. It was obtained using two sublat- 



k^T/J 

B X 



FIG. 2: (a) Heisenberg chain of 8 spins. Nearest-neighbor 
entanglement as a function of magnetic field B and tempera- 
ture T. (b) The same for an Ising spin chain. Here fcs is the 
Boltzmann constant, J and Jx are coupling constants. Light 
color indicates the region where entanglement is detected by 
our method. 



0.77N which increases linearly with the number of spins. 
As shown in Refs. JJ., J/T], when B — the concurrence of 
the two-qubit reduced density matrix in the thermal state 
is obtained as C = max[-((ff//)/iV -I- l)/2,0]. Hence 
C > if (Hh) < ~N and Eh. scp coincides with the 
energy bound for nonzero concurrence. 

Let us now consider the case B > 0. Fig. [3 a) shows 
the nearest-neighbor entanglement vs. B and T. The 
entanglement of formation was computed from the con- 
currence [l5j. Light color indicates the region where the 
thermal ground state is detected as entangled. There are 
regions with C > which are not detected. However, it 
is clear that when the system contains at least a small 
amount of entanglement (~ 0.07) the state is detected as 
entangled. Note that the sharp decrease of the nearest- 
neighbor entanglement around Bcrit = 4 for T = is due 
to a quantum phase transition. 

Another important question is how the temperature 
bound Te depends on the number of particles. For the 
Heisenberg model of even number of spins with B = 
this temperature decreases slowly with N and saturates 
at T « 3.18. Ref. [I3I finds the same bound temperature 
for nonzero concurrence for an infinite system. 



B. XY model 

For the XY Hamiltonian on a d-dimensional cubic lat- 
tice with periodic boundary conditions 

N 

H^y := Y. ^^-i'^-i' + ^.<^4'^ + ^ E -i'^ (10) 

{k,l) k 

the energy of separable states is bounded from below 



(Hxy) > ExY,scp 



-dNM{l + b^/4:) if 6< 2, 
-dNMb if 6 > 2. 



(11) 

Here Jx/y is the nearest-neighbor coupling along the x/y 
direction, B is the magnetic field, M := max(| Ja;|, \Jy\) 
and b := \B\/M. This bound is simply the mean- field 



• — '-'X'^ 



^x ~^'^y^y ^y 



tices and minimizing /xy(s • 
B{sf + sf)/2. ^ 

A one-dimensional spin-1/2 Ising chain is a special case 
of an XY lattice with J^ = 1 and Jy = 0. Fig.j^Jb) shows 
the nearest-neighbor entanglement as a function of B and 
T for this system. According to numerics, Te (computed 
for _B = 1) decreases with increasing N. For N = 00 we 
obtain Tb « 0.41 ^. 



C. Heisenberg coupling between all spin-pairs 

From a theoretical point of view, it is interesting to 
consider a system in which the interactions are described 
by a complete graph rather than a lattice [13 • For the 
following Hamiltonian of N spin-1/2 particles 



Hs 



J^ + J„ + Jj, 



(12) 



the expectation value for separable states is bounded |2C 



Here J, 



x/y/z 



(Hs) > Es.sep ■■= 27V. (13) 

= ^j, a I I and for simplicity N is taken 
to be even. 

Now we could not use the method for partitioning the 
spins into sublattices. The proof of Eq. (|13|l is based on 
the theory of entanglement detection with uncertainty 
relations [23, |211 . For separable states one obtains |2Qj 



(AJ,)2 + (AJ,)2 



(AJ, 



E^ 



E 

k ^ 



(AafOf + (A<0/ + (Aar) 



(fe)N2 



.(fc)^2 



>iV-Ls, 



(14) 



where index I denotes the Z-th subensemble and L5 = 
vai.M{{IS.ax)% + (Acrj,)| -I- (Acr^)|] = 2. Hence Eq. ^ 
follows. The measured energy even gives information on 
the entanglement properties of the system. Based on 
the previous considerations, it can be proved that for 
pure states {Hs) /2 is an upper bound for the number 
of unentangled spins, N^- For mixed states of the form 
P = I]fcPfel*fe)(*fcl we obtain {Hs)/2 > J2kPk^u,k- 
Here N^^k corresponds to the fcth pure subensemble. 

Following the approach of Ref. [l^, the concurrence 
can be computed as a function of the energy. For 
even N the concurrence is C = max{ — [(i/g) + N{N — 
4)]/[2iV(iV-l)],0}. Since for all quantum states (Hs) > 
0, the concurrence is zero for any temperature if A^ > 4. 
Thus our condition can detect multi-qubit entanglement 
even when two-qubit entanglement is not present. 

The thermodynamics of Hs can be obtained by know- 
ing the energy levels and their degeneracies [22| 



E, - 2j-(j + 2), 

(2j-fl)2 / N 



dj = 



A^/2 + ,7 + 1 \N/2 + ]J ' 



(15) 



where < j < N/2. Approximating the binomial in Eq. 
(|15|l by a Gaussian and taking the hmit N ^ oo while 
keeping T/N constant we obtain {Hs) ~ 3NT/{T + 2N) 
and Te ~ 4iV which is in agreement with our numerical 
calculations. Thus Te increases linearly with N, the rea- 
son being that the number of two-body interaction terms 
increases quadratically with the system size. 



value of particle number is constrained to N^. Note that 
for our calculations we assumed that there is at most a 
single atom per lattice site. 



E. Physical realization 



D. Bose-Hubbard model 

Consider now a lattice model, the one-dimensional 
Bose-Hubbard model, in which the number of particles 
can vary on the lattice sites. We use the language of 
second-quantization. Each lattice site corresponds to 
a bosonic mode with a destruction operator a^. The 
Hamiltonian is l24l 



Hf 



{k,i) 



aua'i 



U^alalakttk, 



(16) 



where J is the inter-site tunneling and U is the on-site 
interaction. Let us consider the case when there is at 
most a single particle per lattice site {U ^ J) [23 • Then, 
for separable states the energy is bounded from below as 



(Hb) > Eb, 



-2JNb 1 



1% 

N 



(17) 



where N is the number of lattices sites and iVf, := 
{J2k 'vt'^k) is the number of bosonic particles. For iV == 10 
and Nb = N/2 (half fiUing) we obtain Te ~ 0.69 J. 

Eq. H17|l can be proved as follows. Let us consider 
a site in a pure state |^) = a |0) + (3\1) such that 
|ap -I- |/9p = 1. For this single-site state \{a)\ — \af3\ 
and (atfl) = \(3\^. Hence \{a)\^ = {a''a){l - {a''a)). Now 

using J2{k.i){o-l)M + ^-c- < '^J2k l("fe)P '^^^ can show 
that Eb,scp is an energy bound for product states. It is a 
bound also for mixed separable states of the form Eq. Q 
since EB,scp{Nb) is a convex function. 

Remarkably, the energy minimum for separable states 
equals the minimum for translationally invariant product 
states. In other words, it equals the energy minimum ob- 
tained from the Gutzwiller ansatz |23j if the expectation 



The above methods can be used for entanglement de- 
tection in the following ways: (i) Energy can be di- 
rectly measured in some systems (e.g., optical lattices of 
cold atoms [2j| when used to realize the Bose-Hubbard 
model), (ii) The temperature can be measured and used 
for entanglement detection, (iii) The expectation value of 
the Hamiltonian can be obtained indirectly if the correla- 
tion terms of the Hamiltonian are measured. For exam- 
ple, average correlations X]fe(^a "a )/^; '^ = 2;,y, z 
can be measured in a Heisenberg chain realized with two- 
state bosonic atoms [23. From these correlations (Hh) 
can be computed. 



IV. CONCLUSION 

In summary, we used the Hamiltonian for witnessing 
entanglement in spin models. We also considered bosonic 
lattices. Our further results concerning this system will 
be presented elsewhere [23 • While our method works for 
non-equilibrium systems, we have shown that entangle- 
ment can efficiently be detected by measuring energy in 
a thermal equilibrium. 

Note added. — We presented the idea of using Hamil- 
tonians as entanglement witnesses on a poster during the 
Gordon conference in Ventura, USA in February 2004. 
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